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Abstract 

We explore the effect of varying drive on metastability features exhibited by the vortex matter 
in single crystals of 2H-NbSe2 and CeRu2 with varying degree of random pinning. The metastable 
nature of vortex matter is reflected in the path dependence of the critical current density, which in 
turn is probed in a contact-less way via ac-susceptibility measurements. The sinusoidal ac magnetic 
field applied during ac susceptibility measurements appears to generate a driving force on the vortex 
matter. In a nascent pinned single crystal of 2H-NbSe2, where the peak effect (PE) pertaining to the 
order-disorder phenomenon is a sharp first order like transition, the supercooling feature below the 
peak temperature is easily wiped out by the reorganization caused by the ac driving force. In this 
paper, we elucidate the interplay between the drive and the pinning which can conspire to make the 
path dependent ac-susceptibility response of different metastable vortex states appear identical. An 
optimal balance between the pinning and driving force is needed to view the metastability effects 
in typically weakly pinned specimen of low temperature superconductors. As one uses samples 
with larger pinning in order to differentiate the response of different metastable vortex states, one 
encounters a new phenomena, viz., the second magnetization peak (SMP) anomaly prior to the 
PE. Supercooling/superheating can occur across both the PE and the SMP anomaly and both of 
these are known to display non-linear characteristic as well. Interplay between the path dependence 
in the critical current density and the non-linearity in the electromagnetic response determine the 
metastability effects seen in first and the third harmonic response of the ac susceptibility across the 
temperature regions of the SMP and the PE. The limiting temperature above which metastability 
effects cease can be conveniently located in the third harmonic data, and the observed behavior can 
be rationalized within the Beans Critical State model. A vortex phase diagram showing the different 
vortex phases for a typically weakly pinned specimen has been constructed via the ac susceptibility 
data in a crystal of 2H-NbSe2 which shows the SMP and the PE anomalies. The phase space of 
coexisting weaker and stronger pinned regions has been identified. It can be bifurcated into two 
parts, where the order and disorder dominate, respectively. The former part continuously connects 
to the reentrant disordered vortex phase pertaining to the small bundle pinning regime, where the 
vortices are far apart, interaction effects are weak and the poly crystalline form of flux line lattice 
prevails. 

PACS numbers: 74.25. Qt, 64.70.Dv, 74.25.Dw, 74.25. Sv 
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I. INTRODUCTION 



The richness of the phenomenon of the Peak Effect (PE) in the mixed state of type-II 

nnn 

superconductors has engrossed the vortex physics community for over forty years |1|, I2i |3J, 
3, Q] . The ubiquitous PE phenomenon is an anomalous enhancement in the critical current 
density (J c ) of a superconductor as a function of applied field or temperature in the vicinity 
of the superconducting to the normal state transition The vortex matter can be 

viewed as an elastic medium in a random pinning environment along with the influence 
of thermal fluctuations acting on the pinned vortices. According to a heuristic argument 
due to Pippard 6\ and subsequent theoretical work by Larkin and Ovchinnikov 7], the PE 
is considered to be triggered by an incipient softening of the elastic vortex lattice as the 
normal state is approached. The softer lattice presumably gets conformed to its pinning 
environment more snuggly, thereby producing an enhancement in the pinning (or J c ). The 
notion that the PE can be associated with an inevitable and eventual softening of the 
vortex lattice encouraged the widespread belief that the PE is a (precursor) signature of the 
phenomenon of flux line lattice (FLL) melting 0,0. It is useful to recall at this juncture 
that FLL melting phenomenon has been unambiguously verified only in the anisotropic high 
temperature cuprate superconductors (HTSC) [10]. The PE phenomenon has, however, been 
widely studied in low T c superconductors (LTSC) as well as in the HTSC. In the LTSC, 
the smallness of the Ginzburg number, which measures the relative importance of thermal 
fluctuations vis-a-vis superconducting condensation energy, makes the thermally driven FLL 



melting line lie very close to the H c2 (T) line |3j. The investigations in a variety of samples of 
LTSC reveal that the anomalous behavior pertaining to the PE could surface up sufficiently 
below the H C 2(T)/T C (H) values and the separation between the onset position of the PE 
anomaly and T C (H) line correlates well with the level of effective disorder in the sample 



11 1 . Such observations seem to suggest that, generically, the PE as a phenomenon needs a 
clarification and distinction from the pristine thermally driven FLL melting transition. 

In any realistic sample of a type II superconductor, the inevitable presence of residual 
quenched random disorder is anticipated to compromise the perfect translational symmetry 
of the Abrikosov flux line lattice state. However, it was argued by Larkin and Ovchinnikov 
(LO) 12j that in the presence of pinning, the spatial order in the vortex matter could survive 



within a domain having dimensions much larger than the intervortex spacing, but smaller 
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than the typical sample size. In recent times, Giamarchi and Le Doussal have shown 
that the spatial correlations in the weakly pinned vortex matter could decay much more 
slowly, in an algebraic manner, such that a good positional and high orientational order, 
equivalent to a quasi-long range order (notionally a 'Bragg glass' state), can be observed 
over length scales comparable to the sample dimensions. A contemporary view that has 
gained acceptance about the PE anomaly is that it represents a transition from a weakly 
pinned Bragg glass state to a stronger pinned multi (or micro) domain vortex glass state. 
Bulk ac susceptibility studies in a very weakly (nascent) pinned superconducting specimen 
of 2H-NbSeo had pointed towards the association of the PE with the first order nature of the 

n 

underlying transition in the pinned vortex matter |14j . Eventually, the local ac susceptibility 
measurements by scanning micro-Hall bar microscopy j^J in the typically weakly pinned 
single crystals of 2H-NbSe2 directly elucidated the presence of an interface separating the 
weaker and stronger pinned regions across the peak effect region, thereby endorsing its first 
order nature and the associated metastability effects. 

The presence of a first order transition line in the magnetic phase (H, T) diagram of a 
weakly pinned superconductor imbibes in it the notion of supercooling and superheating 
effects across it y y y. An STM imaging study [l8| alongwith an analysis of the 
instantaneous positions of the individual vortices in a very weakly pinned crystal of 2H- 
NbSe2 led to the surmise that neighboring vortices execute collective motion below the onset 
temperature of the PE, which transforms to positional fluctuations (random excursions) 
of individual vortices above it. A small angle neutron scattering (SANS) study |3] in a 
typically weakly pinned crystal of Nb could reveal the superheating of the ordered Bragg 
glass phase above the onset temperature of the PE alongwith the usual supercooling of the 
disordered vortex glass phase below it. An elucidation of the superheating effects across 
the PE in a variety of single crystals of 2H-NbSe2 by bulk transport measurements using a 
fast current ramp procedure [2(j] has resulted in the demarcation of a generic spinodal line 
corresponding to the limit of the superheating of the ordered state. Above the spinodal line, 
the threshold force needed to depin and drive the lattice (i.e., critical current density J c ) is 
independent of the thermomagnetic history of the underlying pinned vortex matter. 

The work of Xiao et al. |20( on the spinodal line in conjunction with the correlation 
re effective disorder, the structure across the PE and the metastability effects 
raises an issue relating to the possible connection between the quenched random 
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disorder and the details characterizing the spinodal nature. In the present rep ort, we shall 



focus on the results pertaining to this issue. It had been shown earlier |14l . Il7t |2l| that 
with the gradual enhancement in the quenched random pinning, the first order nature of 
the order- disorder transition associated with the PE is not severely compromised; only the 
process of disordering starts to comprise several distinct steps. In a nascent pinned sample 



14j . the interface separating the weaker and stronger pinned regions is perhaps very fragile, 
and the metastability effects in such a sample get wiped out easily by an infinitesimally 
small driving force. In typically weakly pinned specimen of LTSC systems (like, Nb crystal 
used in SANS experiment [XQ| , or the 2H-NbSe2 crystals used in micro-Hall bar microscopy 

^]), the metastable states are more robust and the thermomagnetic history effects can be 
conveniently probed via contactless ac susceptibility measurements, with a minimal level 
of an ac driving force. By performing broad band noise measurements in the in-phase ac 
susceptibility (xL) an< ^ the third harmonic (x&J measurements, one can demarcate the limits 
of metastability and determine the spinodal line in the vortex phase diagram. We elucidate 

;he above stated behavior in the samples of two very different superconducting compounds 

2l| . viz., 2H-NbSe2 (T c ~ 6 K) and CeRu2 (T c ~ 6.3 K) to establish the generic nature of 
the spinodal line, so determined. 

II. EXPERIMENTAL 

The ac susceptibility studies have been carried out using a conventional double coil ar- 
rangement affixed coaxially inside a superconducting magnet coil and employing a mutual 



inductance bridge 22j. Most of the data have been recorded in the frequency interval of 21 



Hz to 211 Hz and with an ac amplitude (h ac ) lying in the range of 0.5 Oe to 3.5 Oejr.m.s.). 
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2811 and a 



The samples studied include the crystals X, Y', Y, Z and Z' of 2H-NbSe2 
single crystal of CeRu2 (21 ] 



]. The crystal X of 2H-NbSe2 is the most weakly pinned sample, 
with a T c (0) = 7.225 K, the crystals Y' and Y are somewhat stronger pinned than the crystal 
X with T c (0) = 7.25 K and 7.17 K, respectively. The crystals Z and Z' have T c (0) values in 
the neighborhood of ~ 6 K and they probably contain few hundred ppm of Fe impurity in 
them |21[. From the sample X to Z', the J r va lues at 4.2 K in the low field region (~ IkOe) 
vary from about 10 A/cm 2 to 10 3 A/cm 2 [23]. We reckon that the FLL correlations in the 
vortex matter in the crystals Z and Z' of 2H-NbSe2 are like those in the typically weakly 
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tinned crystals of other low T c superconducting compounds, e.g., CeRu2, Ca 3 Rh4Sni 3 , etc. 
2l|, 124 ). 



III. RESULTS AND DISCUSSIONS 

A. AC susceptibility measurements in the crystals of 2H-NbSe2 with varying pin- 
ning 

Figure 1 shows the in-phase ac susceptibility response, xL> m * ne different crystals X, Y', 
Y, Z and Z' of 2H-NbSe2 in a field of 15 kOe. The ac susceptibility response varies as [25], 
Xu ~ — 1 + oih ac /J c , where a is a geometry and size dependent factor. This implies that 
the xL response closely mimics the behavior of the J c in the sample. In different samples, 
we compare the J c values of the vortex states corresponding to different thermomagnetic 
histories, viz., the zero field cooled (ZFC) state and the two field cooled (FC) modes. The 
response of a FC state is first measured while cooling down (FCC) , and later measured while 
warming up (FCW). 

Figure 1(a) shows the xL response in the cleanest crystal X for the vortex states in the 
ZFC, FCC and FCW modes. All the three responses are observed to overlap, such that no 
distinction can be made amongst them. In this crystal, we find the PE to be a very sharp 
feature (which occurs over a width of 20-40 mK or so). Considering that the crystal X has 
the least amount of disorder, one may argue that the properties of the first order nature of 
the PE transition should get best exemplified in this crystal. It was, therefore, expected that 
a substantial hysteresis in the responses should have been evident between the ZFC and 
the FCC states due to the possible supercooling effect below the onset temperature of the 
first order transition. While the data shown in figure 1(a) for the vortex states in a field of 15 
kOe is a representative one for the crystal X, we found the absence of any thermomagnetic 
history dependence in the xL response (i.e., differences between the ZFC and the FCC or 
the FCW states) across the PE transition at all magnetic fields in this sample (all data not 
being shown here). A possible reason for the absence of the hysteresis between the ZFC 
and the FC states in the sample X is a subtle effect of the superimposed ac magnetic field 
h ac . We elucidate next an additional facet of h ac in the ac susceptibility measurements on 
the pinned vortex matter, which is different from its obvious role in measuring the shielding 
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response from the macroscopic screening currents induced by its imposition on the sample. 

Figure 1(b) displays the xL response of the vortex state with two different thermomagnetic 
histories, viz., the ZFC and the FCW, in the crystal Y' at 15 kOe. We observe that the 
sharpness of the PE in this crystal is comparable to that in the crystal X, but the response 
of the ZFC and FC vortex states can be made to differ, unlike the situation observed in the 
crystal X. We find that for the ZFC and FCW modes, the xL responses in the crystal Y' 
coincide if the h ac is kept 'on' while the sample is initially cooled down in a dc magnetic field 
to the lowest temperature (which was 4.3 K in the present case). Such a situation is identical 
to that in the sample X. However, if the sample Y' is field cooled with the h ac switched off 
during the cooling, and it is switched back on only when one begins the FCW measurement, 
then one observes a measurable difference in the ZFC and the FCW responses. This behavior 
is indeed different from that observed in the crystal X. These imply that depending on the 
pinning in the sample, h ac has the potential to modify the path dependence in J c . Such 
an attribute becomes further exemplified as we examine the results on the samples with 
higher pinning, viz., the crystals Y and Z'. In these samples, it is significant to note that 
the differences in the responses of the ZFC and the FCW states survive, whether the h ac 
is kept on or off during the initial field cooling procedure (see figures 1(c) and 1(d)). Note 
that the limiting temperature below which the xL response exhibits thermomagnetic history 
dependence (cf. the positions of arrows marked in the different panels of figure 1) could be 
a function of h ac and the extent of effective pinning in the sample. To emphasize another 
observation, the notional peak position of the PE in the ZFC (or the FCW) mode does not 
necessarily mark the limiting temperature above which the thermomagnetic history effects 
in the xL response cease. 

^From the xL responses in figures 1(c) and 1(d), it is evident that yj c < J[ cw . To under- 
stand the behavior of J c corresponding to the vortex states with different thermomagnetic 

n 

histories, we take recourse to a contempory view |4| that in the PE region the vortex state 
comprises an admixed inhomogeneous phase, with the coexistence of ordered and disordered 
regions. The ZFC state prior to the PE represents an ordered weakly pinned state, which is 
characterized by a low J c value. The FCW state on the other hand could correspond to the 
phase above the peak temperature, viz., a predominantly disordered, strongly pinned vortex 
state that can get supercooled down to the low temperatures during field cooling. The FCW 
state, therefore, carries a higher J c value. We have elucidated through the results in figure 
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1(b) that the non-observation of metastable phases in the xL response need not imply the 
absence of path dependence in J c . 

In the presence of very weak pinning, the metastable states of the vortex matter could 
be very fragile. Any attempt to couple the system to the external environment perturbs the 
underlying vortex matter such that it is driven into a different state, thereby masking (i.e., 
altering) the pristine information pertaining to the met ast ability. The process of keeping 
h ac switched 'on' during field cooling the vortex matter can also help the disordered regions 
in exploring the possibilities of transformation into the ordered regions of the ZFC mode. 
This is best exemplified by the data in figure 1(a), wherein in the crystal X, the xL response 
is identical, irrespective of whether the h ac was kept 'on' or 'off during initial field cooling. 
In the sample X, even while h ac was kept off during FC, the application of a small h ac to 
measure the xL response during the FCW mode resulted in driving the field cooled state 
towards the ordered state of the ZFC mode. 

In the sample Y', where the pinning is somewhat stronger than that in sample X, the 
option of keeping the h ac switched on/ off during the initial field cooldown procedure demon- 
strates the competition between progressive enhancement in pinning and the annealing effect 
of h ac . The enhanced pinning in sample Y' arrests the hysteresis across the PE and allows 
the difference between the ZFC and the FCW state to be observed clearly before the onset 
temperature of the peak effect (T° n ). The inevitable dynamical changes which accompany 
the PE phenomenon make the xL responses for the ZFC and the FCW states nearly identical 
for T > T° n in the sample Y'. The progressive increase in the pinning can therefore aid the 
metastability of the supercooled (FC) phase against the thermal fluctuations and the possi- 
ble perturbations from h ac . To search for the true location, where the path dependence in J c 
ceases, in other words to determine the spinodal temperature (T*) which is not influenced 
by the driving forces, one needs to resort to a sample with optimal pinning. Therefore, to 
generate metastability in the vortex matter and to sustain its imprint in a measurement, one 
requires a crucial balance between the strength of the driving force (h ac /transport current), 
disorder /pinning and the temperature. At this juncture, it may also be pertinent to recall 
and compare the results of the electrical transport data of Xiao et al. [26] in a crystal of 
2H-NbSe2 with xL response as in our samples Y and Y'. It has been reported that the 
transport J c value for the FC state in a very weakly pinned specimen differs from that in 
the ZFC state only during the first ramp-up of the current while recording the I-V data. 
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The J c values determined during the ramp-down cycle or during the subsequent ramp-up 
cycles were found to be equal to that for the ZFC state. The passage of transport current 
presumably reorganizes the disordered state of the FC mode towards the ordered ZFC mode, 
when the pinning effects are in the nascent stage. Xiao et al. 20] had to later adopt the 
procedure of fast current ramp to explore the metastable states above the notional onset 
temperature of the PE. They found that the superheating of the ordered state could be 
observed upto a limiting temperature T*, which exceeds the notional peak temperature of 
the PE for the vortex state prepared in the ZFC mode j^j]. 

The plots in figure 1(d) appear to imply that the metastability effects in the yL response 
would cease above the peak temperature T p , as was conventionally believed QHQ The 
FC state in this sample can be ascribed to a freezing in of the disordered vortex state at 
the pea k position of the PE, where the lattice was believed to be homogenously amorphous 
I21I However, the recent work by Xiao et al. (2^] desires a serious revision of the 

above presumption. A transformation that sets in at T° n results in an admixed state of 
ordered and disordered phases. Both these fractions carry different J c values and possess 
different temperature dependences. As the temperature increases from T° n towards T p , the 
fraction of the ordered vortex phase decreases and that of the disordered phase increases, 
correspondingly the observed J c value also increases. The relative fractions of these two 
phases and their temperature dependences determine the peak position (i.e., T p /H p ) of 
the PE anomaly. In view of this, the relationship, if any, between T p and the limiting 
temperature at which the metastability effects cease is not very apparent. We have set out 
to clarify this issue via the ac susceptibility measurements. 



B. Comparison of AC susceptibility data in the crystals of 2H-NbSe2 and CeRu2 

As stated earlier, in order to study the features associated with a first order transition one 
should in principle attempt to select a sample with the least amount of disorder. However, 
the data in figure 1 show that to study features related to the spinodal temperature T* via 
the ac susceptibility measurements, we were compelled to choose the samples of Z and Z' 
category, which have significantly higher level of disorder than the sample X. While it is 
clear that with an enhanced pinning, the metastability features of the vortex matter stand 
preserved, it is not however apparent as to whether the first order character of the PE 
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transition stands retained at comparatively higher levels of disorder. A trend which emerges 
from the xL response in figures 1(a) to 1(d), is that between T° n and T p , the PE phenomenon 
undergoes a broadening with progressive increase in the effective pinning. It could be argued 
that this indicates that a sharp first order like PE transition is being transformed into a 
continuous transition. However, in this context, it is useful to recall that a novel notion 
of ( m ul t i- St ep> fracture had been de_ed ear.,, fl Q. The fr* o rder nature of 
the PE transition is not destroyed in samples with higher pinning, it only transforms into 
multiple small first order like jumps, which comprise the notion of fracturing. 

Figure 2 shows a comparison of the ac susceptibility response at different ac amplitudes in 
a crystal Z of 2H-NbSe2 with that in a crystal of CeRu2 (cubic system with T c (0) ~ 6.3 K), 



with a level of pinning comparable to that in the crystal Z [2l|. In both these samples, we 
observe in figures 2(a) and 2(c) that although the PE is a broad feature, yet it is characterized 
by multiple sharp jumps commencing at T°" in the response, characterizing the 

fracturing transition. A well ordered vortex phase is considered to transform into a multi- 
domain like state, with pockets of ordered and disordered regions created, possibly via the 
progressive permeation of topological defects like dislocations, into the ordered vortex phase. 
At each first order like jump, some particular pockets with a collection of well ordered vortices 
(with a pocket characterized by a low J c value), transform into a disordered phase with a 
higher value of J c , locally. More the number of such domains, with intervening disordered 
vortex phase, the average J c of the sample increases leading to a PE like feature above T°™. 
The coexisting ordered and disordered phases are not nucleated by merely supercooling the 
vortex matter below a first order transition, they can also get generated due to the staggered 
nature of the first order transition as a consequence of the fracturing phenomenon. The 
fraction of domains with disordered vortices is larger while field cooling below the fracturing 
transition, and this result in prominence of the metastability effect. 

The above description finds an echo in the plots contained in the figures 2(b) and 2(d). 
We observe that in the samples of 2H-NbSe2 and CeRu2, there is a substantial difference 
between the ZFC and the FCW states for h ac = 0.5 Oe (cf. figures 2(a) and 2(c)). For 
larger h ac values (i.e., between 2 Oe and 3.5 Oe, as shown in the plots for figures 2(b) and 
2(d), the PE becomes a narrower transition and the fracturing features can get suppressed 
and, consequently the metastability response (as reflected by the differences between the 
ZFC and FCW responses) also stands reduced. The connection between the weakening 
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of the fracturing transition and the decrease in metastability is presumably related to the 
annealing effect of the h ac drive in modifying the relative fractions of the ordered and disorder 
regions between T° n and T p . With a larger h ac drive used for x' u measurements, the residual 
disordered pockets in the ordered vortex matter get transformed into the ordered ones, prior 
to the onset of the PE. The vortex matter thus approaches a single domain like picture of 
the ordered vortices as in the samples X and Y' of 2H-NbSe2, which then undergoes a sharp 
PE transition, rather than a multi-step fracturing transition. It is also pertinent to note 
that with the enhancement in the ordered fraction of vortices in presence of a driving force, 
the disordered domains of vortices can permeate and survive in the midst of the ordered 
vortex matter only at relatively higher temperature (cf. the T° n and T p values in the figures 
2(b) and 2(d) with those in the figures 2(a) and 2(c), respectively). It is also instructive 
to compare FCW responses in the figures 2(a) and 2(c) with those in the figures 2(b) and 
2(d), respectively. The supercooled state in the FCW mode need not undergo order-disorder 
transition across the PE region; such a behavior is evident when h ac is small (cf. figures 2(a) 
and 2(c)). However, a larger h ac drive can reorganize the initial FCC state during warm-up 
and eventually depict an order-disorder transition across the temperature region of PE, as 
evident in figures 2(b) and 2(d). 

Having explained the metastability features associated with the fracturing transition seen 
in samples with higher pinning, and determining the optimum range of parameters for h ac , 
disorder and temperature, we now dwell on the demarcation of the limit of metastability in 
the vortex phase diagram. Figure 3 shows the response at different applied fields for 

ZFC, FCC and FCW modes in the sample Z' of 2H-NbSe2- Note that at the lowest temper- 
ature, the FCW and FCC responses are identical, at higher temperatures the FCW state 
is more ordered in comparison to the FCC state. Apart from the behavior of metastability 
and the PE, we observe an additional anomaly in the behavior of the xL response on the 
ZFC branch. In samples having a level of pinning as in the crystal Z', we can notice the 
occurrence of an anomalous enhancement in J c deep in the mixed state beginning at T°^ p , 
corresponding to the second magnetization peak (SMP) feature observed in the isothermal 
magnetization hysteresis (M-H) loops well before the onset of the usual PE at T° n (2^ . The 
SMP feature is different from the fracturing phenomenon which gets triggered at (or near) 
T° n (see figure 3(a)). It suffices to recall and state here that the (two) anomalies of the 
SMP and the PE are distinct and different j^J. It may be worthwhile now to ask as to at 
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what stage does the phase coexistence of ordered and disordered phase commence, in the 
presence of a SMP anomaly. A conventional notion as alluded to in the description earlier 
would imply that the phase coexistence region commenced from the onset position of the 
PE. We shall now show that such a notion desires a revision in view of the data presented 
in the panels of figure 3. Closely associated with the notion of phase coexistence and the 
first order nature of the PE is that of the metastability. 

While examining the metastability response, it is instructive to focus on the interesting 
behavior illustrated in the data recorded at 6 kOe (cf. figure 3(b)). Note that the peak 
temperature of the PE is different for different thermomagnetic histories, the highest value 
being in the ZFC state, T Z J C . Similar trend can be noted at other magnetic fields (cf. 
figures 3(c) and 3(d)). Such a behavior for the xL response for the vortex states with 
different thermomagnetic histories, is analogous to that reported for the transport critical 
current I C (T) data in a crystal of 2H-NbSe 2 by Xiao et al. (cf. figure 1 of Ref.|20j|). Vortex 
states with different histories are characterized by the different fractions of the ordered and 
the disordered phases. The observation that T^ cc < T^ cw < T^ c implies that the limiting 
temperature at which the sample is homogenously filled with the disordered phase does not 
lie in the temperature interval, from T^ cc to T Z J C . 

The xL(T) response for a given H is essentially dictated by J C (T). The temperature 
above which xL(T) (or, J C (T)) becomes independent of the thermomagnetic history of the 
specimen appears to lie even above the (highest) peak temperature, i.e., at T > T**°. 
Figure 4(a) focuses attention onto the plots showing (x'u ~ X^ w ) an d (x'^ c ~ x'Ji cc ) & t 
12 kOe in the sample Z/. In this panel one can identify the limiting temperature T* above 
which the xL response becomes path independent, viz., where the differences ~ x'w) 

and c — x'Jj cc ) vanish. In the sample Z/, due to the relatively strong pinning and the 
fracturing phenomenon, the determination of the limit of metastability effects is insensitive 
to perturbations from h ac . We believe that the above procedure yields a reasonable estimate 
of the spinodal temperature, T*, for 2H-NbSe2 samples with J c (~ 1 kOe, 4.2 K) > 1000 
Amps/cm 2 . However, the above procedure may not suffice for all weakly pinned samples, 
like, the crystals X and Y' of 2H-NbSe2- It is therefore desirable to explore an alternative 
way to determine the value of T*. 
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C. Spinodal temperature T* and the third harmonic ac susceptibility measure- 
ments 



The critical state model (CSM) relates the macroscopic J c to the hysteretic mag netization 



response of a superconductor 
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301 131| . It is well documented 
features get added to the pristine relationship between the J C (H) and the magnetization 
hysteresis, when J C (H) does not remain single valued function of H. For instance, the minor 
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hysteresis loops can display anomalous characteristics, like, an asymmetric shape 
excursions beyond the envelope M-H loop |3Jj|, etc., when the J C (H) turns path dependent. 
In the back drop of these observations, it is instructive to examine the response of the 
third harmonic of the ac susceptibility across the SMP and PE regions, i.e., from below the 
onset temperature of the first of the two anomalous variations in J c upto the irreversibility 
temperature (T irr ), where the (bulk) J c ceases. To be specific, consider the cool down of a 
weakly pinned type-II superconducting sample from above a given T C (H) value to below its 
Tj rr (H), where the finiteness of J c would result in a non-linear magnetization response which 
could generate a measurable third harmonic signal in the ac susceptibility measurements. 
Such a third harmonic signal would be expected to follow the increase in J C (T) for a given 
I as (Tirr - T) increases, as per a prescription of the CSM for path independent J C (H, T) 
30|,|31 1 . The onset of the history dependence in J C (H) could compromise the above stated 
notion, arising from the applicability of the CSM. Below T*, due to the onset of the path 
dependence in J c , qualitative changes could occur in the behavior of the nonlinear response 
at the nucleation of an inhomogeneous phase. 

The plots in figures 4(b) to 4(f) explore the limit of the path dependence in J C (H) via the 
study of the temperature dependence of the third harmonic ac susceptibility data, x 3[J , i n the 
sample Z' of 2H-NbSe2. In figure 4(b), we draw attention specifically to the behavior of 
just above T*. Warming up from the low temperature side, as the temperature crosses the 
limit T*, xkj(T) is seen to monotonically decrease and eventually vanish at the irreversibility 
temperature, T irr (< T C (H)). There does not appear any simple correspondence between the 
XU T ) and xL( T ) for T™„ < T < T* (cf. figures 3(d) and 4(b)). Figure 3(d) shows that 
the xL(T) monotonically decreases above T^ c , reflecting the collapse in J c above the peak 
position of the PE. The X 3w (T), on the other hand in figure 4(b), appears to enhance between 
T Z J C and T*. The xk;(T) response is seen to turn around only above T*, and follow the J C (T) 
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thereafter. Thus, the onset of an inhomogeneous phase at T < T*, produces modulations in 
the behavior of nonlinearity or in the x&j response. 

To establish the assertion on the imprint of the limit of the path independence in J c in 
the X3^(T) data, we can examine the above stated behavior at different fields in figures 4(c) 
to 4(f), and note one to one correlation between values of T* determined from xLO-0 data 
in different thermomagnetic histories and the limiting temperatures above which xk,(T) 
monotonically decreases. The T° n , T^ c and T* values have been marked for each of the 
curves in figure 4. At T < T*, the modulations in X3aj(T) display complex behavior, dictated 
by the (intrinsic) changes in the relative fraction of the ordered and disordered phases (path 
dependence in J C (T)) and the additional changes induced by the h ac value in which the 
Xg w (T) data are recorded. The T*, however, does not appear to vary in any noticeable 
manner with the amplitude of h ac (all data not shown here). It is also worthwhile to note 
that the modulation in the xL,(T) signal ceases close to T°^ p . In fact, if we approach T°^ p 
from the lower temperature side, the xk^T) plots appear to show an enhancement in the 
non-linear response at T > T°^ p . Such an observation would get rationalized by invoking 
the notion that the SMP and the PE are non-linear phenomena In the absence of 
complications induced by (possible) phase-coexistence and history effects, the xkj(T) would 
have been expected to follow the modulation as displayed by xL(T)- We believe that the 
possibility of phase-coexistence commences at T°™ and lasts upto T*. 



D. Correlation between X3^(T) and measurement of noise in xL(T) 

To further substantiate that the values of T* determined from the xL, response indeed 
physically correspond to the spinodal temperatures, we explore the correlation between the 
X 3a; (T) and the noise signal in xL(T) 13 , which can be easily recorded using a Lock-in am- 
plifier having a wide band filter option. The above stated noise signal is believed to measure 



the fluctuations in xL(T) and it is argued |2lJ to reflect the possibility of transformations 



amongst metastable states accessible from a given mode (ZFC/FC). 

Figure 5 depicts the plots of X 3 2 J C (T) (in h ac = 3.5 Oe (r.m.s.)) and the noise in x^ c (T) 
in h ac = 0.5 Oe (r.m.s.) in a field of 5 kOe in the sample Z of 2H-NbSe2 and in field of 13.5 
kOe in the crystal of CeRu2- In different panels of figure 5, we have marked the respective 
values of T* as determined from the onset of nonmonotonic modulations in x'zw response. 
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In view of the notion of the existence of a homogenously disordered phase above T*, it is 
indeed not a coincidence that the noise signals in the panels (b) and (d) of figure 5 recede 
to the background value as T — *• T*. Taking cue from earlier studies of noise in transport 



experiments Banerjee et al. 



21 1 had argued that the increase in noise at T = T c 



(equivalent to the temperature of the onset of plastic flow, T p i, in Ref. 21 1) reflects the 
possibility of enhancement in transformations amongst the coexisting ^| metastable states 
in a fractured (partially disordered) vortex solid (which exists at T < T*). The setting 
in of the decrease in the noise signal as T — > T p was considered to imply the effect of 
phase cancellation of a large number of incoherent fluctuations as the vortex matter moves 
towards the fully disordered state. In such a framework, the possibility of coexistence of 
ordered pockets embedded in the disordered medium would cease as T — > T*, and the noise 
signal would be expected to recede to the background value. In the context of our present 
results, the T*(H) values indeed represent the notion of the spinodal line, T S (H) 



E. Magnetic Phase Diagram in a typically weakly pinned crystal of 2H-NbSe2 

To comprehend a variety of results presented above, we have constructed a field- 
temperature (H, T) diagram for the sample Z' of 2H-NbSe2 in figure 6. It includes the 
values of the onset temperatures of the SMP and the PE anomalies, along with the values 
of T^ c , T*(H), T C (H) and H p ; aiecm (T). In sample Z' (or Z), there is adequate perceptible 
difference in T Z J C and T*(H), and this difference enhances as H increases. However, in the 
samples Y and Y' of 2H-NbSe2, where the PE manifests as a sharp transition, T* values lie 
in very close proximity of the corresponding peak temperature values. In the vortex phase 
diagram of figure 6, the plateau line, H piate(m (T), passes through the limiting fields Il7| 
below which the collective pinning regime gives way to the small bundle pinning regime. It 
seems appropriate to identify the (H, T) space between the H piateau (T) and the T°" p (H) as 
the Bragg (elastic) glass region 13|. Above the T*(H) line, where the metastability effects 
cease, the vortex matter exists in the pinned amorphous phase. The phase diagram demar- 
cates the regime of phase-coexistence j^fl]. We propose that T°^ p (H) line corresponds to 
the boundary where pockets of disordered phase start proliferating in the ordered vortex 
state, however, the fraction of the ordered phase exceeds that of the disordered phase upto 
T° n . Above T°™, the balance starts to rapidly tilt in favor of the disordered regions and the 
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strongly pinned phase starts determining the overall electromagnetic response. Above the 
spinodal line T*(H), the phase coexistence and metastability features cease and the sample 
is homogenously filled with the disordered regions. The bifurcation of the coexistence phase 
space into regions I and II, where the order and the disorder dominate, respectively, finds 
additional support from the temperature dependence of X^(T) in the FCW mode in figures 
3(c) and (d). Both sets of data show that J C (T) values for the stronger pinned vortex solid 
continue to decrease between T°^ p and T° n . It is only above T° n that the J C (T) values for 
ZFC and FCW modes start to increase. 



IV. SUMMARY 

To summarize, we have investigated the effect of drive and disorder on supercool- 
ing/superheating effects across the anomalous variations in critical current density in several 
single crystals of 2H-NbSe2 and compared their results with those in a single crystal of CeRu2. 
We have mapped out in typically weakly pinned crystals of 2H-NbSe2 and CeRu2 the phase 
coexistence regime of the stronger and the weaker pinned pockets, where metastability ef- 
fects manifest in a prominent manner. The termination of the phase coexistence regime at 
higher field and higher temperature end can be conveniently located in the third harmonic ac 
susceptibility data. At lower fields (< 1 kOe), the region I of the phase coexistence regime, 
in which the order dominates, can be seen to continuously connect to the so called reentrant 
disordered phase 11, [l2|, where the intervortex spacing presumably exceeds the range of 



interaction (ao > A, where A is the penetration depth). Bitter decoration [41] and simulation 
studies Q shown that the reentrant disordered phase comprises polycrystalline form 
of flux line lattice. So long as the domains are large, with radial correlation length R c ^> 
ao, the weaker pinned regions dominate. The electromagnetic response of the sample and 
the imposition of the external driving forces can shrink the stronger pinned domain wall 
regions. However, as the domain sizes shrink on crossing over to the region II of the phase 
coexistence regime, the external driving forces aid the process of complete amorphization of 
the vortex matter. The vortices remain pinned in the amorphous region between T* and T c 
values in the typically weakly pinned crystals of low temperature superconductors. 
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FIG. 1: In-phase ac susceptibility data showing peak effect (PE) in a field of 15 kOe with different 
thermomagnetic histories in single crystals of 2H-NbSe2- The crystals X, Y', Y and Z' have 
progressively enhanced pinning. The crystal Z' has a lower T c (0) value of ~ 6K as it presumably 
contains 200 ppm of Fe impurity. The panel (a) shows that in the most weakly pinned crystal 
X, all the four responses, ZFC, FCC, FCW (recorded after h ac on or off during field cooling) are 
identical. In panel (b), the two FCW responses (with h ac on and off during field cooling) in the 
crystal Y' are different. While the former coincides with ZFC (dotted curve), the latter overlaps 
with FCC (the solid line). The panel (c) shows that in the crystal Y, the xL response in the FCW 
with h ac on differs from that in the ZFC state20In the panel (d), the xL data in the crystal Z' is 
such that the FCW response represents the supercooling of the vortex state existing at the peak 
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FIG. 2: Comparison of the xL data at the fields indicated in the ZFC and the FCW modes with 
different amplitudes of the ac driving force (h ac ) in a typically weakly pinned single crystal Z of 
2H-NbSe 2 (T c (0) « 6.0K) and in a crystal of CeRu 2 (T c (0) « 6.3K). The positions of the onset 
and peak temperatures have been identified in each of the panels. 



21 



4.4 



4.8 



5.2 



5.6 



0, 
-0 
-0 



"D 0. 

CD 
(J) 

03 

E-o. 



o - 



8^ 



-0. 

0. 



8 -: 
o - 



-0.8 



— I — I — I — I — I — I — I — I — I — I — I — I — I — I— 

2H-NbSe 2 ( # Z ' ) ^ on( c 



'(a). 




1 p f 2.6 kOe" 

' ' L 



-I 1 1 1 1 1 1 1 i | 1 r- 

' ,,,„> T on (zfc) 
T on(zfc) J 

ZFC smp t 




FCW^^ 
( b ) FCC 4 



i i i i i i i i i 



"i — 1 — 1 — 1 — 1 — i — r 

-on (zfc) 



5 6 kOe 

1 ' ' I ' 



g " ZFC CTsMpJpL PE 



I 1 1 



FCW 



(c) 



FCC 



H — I — 1- 



| I I — I — h 

-on (zfc) 



ZFC smp| SMP ♦ p 



FCW 
- (, g I FCC 




3.0 3.5 4.0 

T(K) 



4.5 



FIG. 3: xL d a * a at different applied fields for different thermomagnetic histories in a crystal Z' 
(T c (0) ~ 6.0K) of 2H-NbSe2- Note the occurrence of two anomalous variations corresponding to 
the SMP and the PE in each of the panels. The xL(^) curves at H = 2.6 kOe in the panel (a) 
display the notion of stepwise fracturing across the PE, as in the Fig. 2 (a) for the sample Z at 5 
kOe. The onset temperature of the SMP and the PE and the peak position of the PE in the ZFC 
mode have been marked in all the panels. 



22 



3.0 



3.5 



4.0 



4.5 



< 



0.4 - 



0.0 



Tzfc 
T* 

I V - V ) >X iP 1 

( a ) 2H-NbSe, ( # Z ) ^* 



0.0 - 
-0.1 - 



i 1 1 1 1 1 1 1 1 1 r- 

-p on * i zfc i fee 




12 kC^ 



_i i i i I i i i i_ 



(b) 



"I 1 1 1 1 1 1 1 T" 



_l I I l_ 



i 1 1 r~ 

>T zfc 



-p on 
| smp 



H || C T °n / \'l C 

I . I . . T. ' 



"S 0.00 - 
N 

2 -0.08 - 

£ o.o - 

o 
c 



-1 1 1- 



(C) jon 

2H-NbSe„ ( # Z' ) ♦ p 



'TV 



|T* 3.0kOeH 



(d) 



\ — i — | — i — i — i — i — i — i i J . i — i — i — 

1 -r on .■ * 



3.6 kOe 
H || c 1 




n ,9.0 kOe 2H-NbSe 2 (#Z') 

H 1 1 1 1 1 1 1 1 — I 1 1 1 1 1 h 




FIG. 4: The panels (a) and (b) show the plots of difference susceptibility i.e., ( x'u - 

Xui ^ CC ) and ( x'w - ^ CW )i an d the third harmonic signal x'zui m a °f 12 kOe (H\\c) in the 
crystal Z' of 2H-NbSe 2 . The positions of T s °" p , T° n and T p for the ZFC mode, as evident from 
the curve in Fig. 3 (d), have been identified in these two panels. T* in the panel (a) identifies the 
temperature at which the difference (xL - Xu ^ CC ) vanishes and merges into the baseline. This 
value of T* is then marked in the panel (b). Note that x'zuj res P° nse shows a monotonic decrease at 
T > T*. The panels (c) to (f) show x%u> curves at different applied fields with h ac = 2.50e (r.m.s.) 
in the crystal Z' of 2H-NbSe2- The respective T* values have been identified in these panels. It is 
apparent that the x'^w collapses above T* and Merges into the baseline as T — > Tj rr (< T c ), where 
T;rr is the notional irreversibility temperature for a given hac and the frequency value in which 



0.0 
M -0.2 
~ -0.4 



4.8 5.0 5.2 

5kOe 
h = 3.5 Oe ( r.m.s. ) 

'•^—t iT on<zfc) i 
K' T 

■ P ■ 

■ (a) \ iLy 

2H-NbSe, > / T* 

(#z) * 



/t 



* 20 ^ 

to 
•*-> 

o 

\ 10 ^ 



V 

tn 



2H-NbSe 2 

(#Z) ; Tpi l 

■ ■ ■ ■ ■ 

(b) ' WT " 

, ] Br*? 

V^W 5 kOe 

h = 0.5 Oe ( r.m.s. ) 



4.8 5.0 



13.5 kOe I 
h = 3.5 Oe ( r.m.s. ) 



5.2 



(c) 

CeRiL 



0.00 



-0.04 



-0.08 



CeRu 



13.5 kOe 
2 h = 0.5 Oe ( r.m.s. ) 

'fo- 



rt 



(d) l!i J 



20 



10 



4.8 5.0 5.2 4.8 5.0 

T(K) 



5.2 



FIG. 5: A comparison of third harmonic, x'zuj response and the noise signal in the ac susceptibility 
at the fields indicated in the typically weakly pinned single crystals of 2H-NbSe2 ($Z) and CeRu 2 . 
The X3uj data are recorded in h ac of 3.5 Oe (r.m.s.), whereas the measurements on the noise in xL 
[see Ref. 21] were made with h ac = 0.5 Oe (r.m.s.). The j^ 2 -^ value in the panel (a) corresponds 
to the said value identified in the Fig. 2 (b) for the same h ac of 3.5 Oe (r.m.s.). The T p i value in 
the panel (b) identifies the onset temperature of the PE in h ac of 0.5 Oe (r.m.s.) in the sample Z, 
following the nomenclature as in Ref. 21. T p i{H) values denote the boundary of elastic to plastic 
flow for the driven vortex matter. A comparison of data in panels (b) and (d) show that T* values 
identify limiting temperatures above which x'su> responses monotonically decrease and the noise 
signals in xL recede to the respective baselines. 
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FIG. 6: Vortex phase diagram for H\\c in a typically weakly pinned single crystal of 2H-NbSe-2. 
The diagram comprises the values of onset temperatures of the second magnetization peak anomaly 
(T™ p ) and the peak effect (T° n ), the value of the peak temperature of the PE (T p 2/c ) in the ZFC 
mode, the limiting (spinodal) temperature T*, the T C {H) and the values of the crossover field, 
H p iateau{T)i determined from the normalized plots of critical current density versus normalized 
fields as in Ref. 21. The (H, T) phase space above H p i ateau (T) and below (H) identifies the 
ordered Bragg (elastic) glass region. The intervening phase space between T™ p (H) and T*(H) 
comprises the co-existence regime of ordered (weaker pinned) and disordered regions. T p n (H) 
values sub-divide the co-existence phase space into regions I and II in which ordered and disordered 
pockets dominate, respectively. The region I appear continuously connected to the phase space 
of the so called reentrant disordered regime, in which vortices are in the small bundle collective 
pinning regime. The phase space intervening between T* and T c identifies the pinned amorphous 
region. 
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